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ABSTRACT 


A laboratory study was conducted to examine the effect of 
viscous and capillary forces on breakthrough oil recovery in an 
unconsolidated Ottawa silica sand pack. Displacement tests were 
conducted using distilled water to displace iso-octane and Pembina 
crude. These test results were scaled using the ratio of viscous to 


capillary forces (L,)- 


Displacement tests conducted on iso-octane ( Wo / My ple) 
resulted in recoveries at breakthrough of around 72 percent oil. 
Pembina crude ( Ho / Hy > 1 ) had recoveries at breakthrough about 
43 percent for the same core and conditions used for iso-octane. 
Initial water saturations for both systems were not the same and 


decreased with a change from iso-octane to Pembina crude as the test 


fluid. 


Over the range of the universal rate function studied, it was 
observed that no viscous fingering occured in the floods conducted on 
iso-octane. However, as the rate function increased in value there 
was a decrease in the breakthrough recovery of Pembina crude which was 
attributed to the formation of water tongues at the flood front 


(viscous fingering). 


During the tests carried out with Pembina crude it appeared that 
the wettability changed from a water-wet to an oil-wet system. This 
change was assumed to be a function of the effect of oil composition 


on the solid-liquid surface and was time dependent. 
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INTRODUCTION 


Controlling the displacement of oil by water to achieve the 
maximum efficient and economic recovery from reservoir sands is a 
problem of major importance in the oil industry today. In an attempt 
to explain the increase in recovery observed from the waterflood 
process, teen cited a number of parameters which influence the dis- 
placement phenomena; these being capillary forces, oil viscosity, 
interfacial tension, wettability, flow rate and the presence of gas. 
Dimensional analysis applied to these variables results in normalized 
scaling functions which may be used to relate the laboratory model 


to the actual reservoir or prototype. 


Figure 1 illustrates the theoretical breakthrough recovery as 
influenced by the interplay of capillary and viscous forces in a 
linear horizontal model. Region A to B depicts an increase in oil 
recovery with an increase in a dimensionless rate function. This 
increase has been attributed to a decrease in the influence of 
capillary forces on the recovery process. In the laboratory these 
forces are observed as the inlet and outlet end effects caused by a 
discontinuity of the capillary forces across the core ends. As a 
result of the minimization of these forces at an increased value of 
the rate function, a region of stabilized recovery is attained (B to C). 
The phenomenon of bypassing the oil by water, leaving behind isolated 
pockets of hydrocarbons is termed viscous fingering? (Figure 2). At 
point C (Figure 1) the onset of fingering within the porous medium 
occurs and beyond point D viscous fingering has no additional effect 


on recovery. 


If the onset of fingering could be established, and thereby 
avoided, then maximum recovery could be obtained by the waterflood. 
The objective of this study was to establish the magnitude of a 
scaling factor which might be used to predict the onset of finger- 


ing within a porous medium. 
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LITERATURE REVIEW OF 
THE PHENOMENON OF VISCOUS FINGERING 


Secondary recovery of oil by displacement with water is not a 
simple straight forward process. It involves a complex displacement 
phenomenon which is affected by the reservoir rock and fluid properties. 
The interaction between the capillary and viscous forces is one of the 
phenomena that must be fully understood and controlled for the attain- 


ment of maximum oil recovery. 


The influence exerted by these forces on recovery was illustrat- 
ed in Figure 1. Jones Parra et ike and Kyte et a observed that the 
lower recoveries in the capillary region were due to water fingering 
at the inlet end of the core, during slow displacement rates. As the 
waterflood progressed these fingers coalesced into a sharp planar 
interface. It was, therefore, concluded that a capillary discontinuity 
across the core end existed, known as the inlet end effect, and it was 
this effect, not the viscous forces, which caused the finger formation. 
At the outflow end of the sand column there exists a capillary dis- 
continuity. The buildup of a water saturation in the boundary grain 
layers as a result of the capillary discontinuity decreases the per- 
meability to oil. This laboratory phenomenon is known as the outket 
end effect. This outlet end effect decreases as the scaling factor 


increases in magnitude. 


Where the driving fluid is more mobile than the driven fluid 
the interface between the two fluids may become unstable and lead to 
gross channeling causing a severe decrease in recovery efficiency. It 
was this observation which prompted Engelberts and Klinkenberg to 
call this phenomenon viscous fingering. In all cases where viscous 
fingering has been visually observed a marked decrease in recovery 


was also noted. 


A quantitative analysis of fingering was presented by Chouke” 


in 1958. His mathematical method predicts that fingering, due to 
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viscous forces, will occur if the finger thickness is slightly less 
than the core diameter. By varying the injection rate, oil-water vis- 
cosity ratio and wetting conditions, densa, experimentally substant- 
iated Choukes theory by showing a marked decrease in recovery at the 
predicted value for the onset of fingering. By calculating the onset of 
fingering for the earlier work of Engelberts and Klinkenberg the same 


observation was possible. 


The displacement process is also influenced by the preference 
of a solid to be wetted by a specific liquid, called the wettabihity 
of that material. Pigeon: describes a number of saturation states giving 
rise to various phases of wettability within the two general categor- 


ies of water-wet and oil-wet systems. 


A porous medium is said to be water-wet if the contact angle, 
of the water with the oil on the rock surface, as measured through the 
water phase, is less than 90 degrees and oil-wet if it is greater than 
90 degrees. Since wettability is often dependent upon which fluid first 
contacts the grain surfaces, the majority of reservoirs are water-wet. 

; "However, prolonged contact with oil, 

particularly if the latter contains 

easily absorbable polar compounds 

may partially and locally change the 

wetting properties of a reservoir 

and give rise to wettability 

praddentet) 
An oil-wet sand is more likely to have less connate water saturation 
than a water-wet sand due to expulsion of the water as the oil prefer- 
entially wets the rock surface. Experimental evidence of deHaan 
suggests that fingering occurs primarily in water-wet systems with 
no evidence of water channeling observed in the oil-wet case. 


In order to obtain these macroscopic deformations at the stable 


* 
Pirson, S.J.: "Oil Reservoir Engineering,'' McGraw Hill, 


Toronto, (1958), p.69. 
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interface, force must be applied along the interface. This force acts 
to overcome the interfacial tension which may be defined as the surface 
energy between a liquid and a liquid or a liquid and a solid. For the 
water-wet system this force is much less than for the oil-wet system 
since less energy is required to overcome the interfacial tension be- 
tween the oil and water than between the oil and rock. This implies 
that if viscous fingering occurs in the oil-wet medium it would only 
be manifested at larger throughput rates than in the water-wet case. 
That is the degree of wettability will affect the viscous fingering 
phenomena, and such a phenomenon could occur in an oil-wet system 


contrary to the observations of deHaan. 
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THEORY OF 
IMMISCIBLE OIL-WATER DISPLACEMENT 

Displacement of a wetting by a non<wetting fluid, or vice 
versa, from a linear porous media has been mathematically modelled 
by a number of. researchers. They assume the primary displacement 
mechanism to be the movement of a steep saturation front through the 
reservoir followed by an additional gradual oil displacement known as 
the subordinate phase. These theoretical models have provided a means 
of predicting the performance behaviour of waterfloods within the 


limits of the mathematical assumptions. 


In 1941 Deverere’ evaluated Darcy's law within a gravitational 
field for two liquids, water and oil. By combining these relation- 


ships with the formula for capillary pressure: 


when Pa j=ne - P (1) 
c Oo w 
Re = capillary pressure (atm) 
a = pressure in the oil phase (atm) 
and ie = pressure in the water phase (atm) 


the fractional flow equation was obtained: 


Ko oP. - g Ao sina | 
feta ae i 
qu du 
f= ee Oe eee (2) 
Ag ee 
who 
f = fraction of water in the flowing 
when w 
stream at any point 
Ky = effective permeability to oil (darcies) 
and Kk = effective permeability to water (darcies) 
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eee + q _ = total fluid flow rate 
per unit cross sectional 


area (a apefcne) 


U_ = viscosity of oil (cp) 


L.. = viscosity of water (cp) 


u = flow direction (cm) 


g = gravitational constant = 1033 


(cm of water/atm) 


Ap = density difference (gm/cc) 


and Qa = angle between u and the horizontal 

Later Buckley and hevereres: (1942) undertook a theoretical 
analysis for the displacement of two immiscible fluids. They assumed 
the gravitational and capillary effects negligible and combined the 
fractional flow equation with a material balance to yield the frontal 


advance rate equation: 


Q. of 
du = — | —“|dt (3) 
oA | as. 
when Q, = total displacing fluid entering 


the system (cc/sec) 
@ = porosity (fraction of pore volume) 
2 
A = cross sectional area (cm ) 


S = saturation to water (fraction of 


pore volume) 


and t = time (sec) 


This equation relates the linear advance of an interface of 
constant saturation with an increment of time under the application of 


a constant fluid injection rate through the porous medium. 


Lavaes tan aah cto 
t: 

A t Bie i 
Lspliayasnd & Asoduabee (os0ny 2388 
beinuRad voit ; BbigEs efijee . 


ay ‘beatdmog bie» oidtaeh si —_ 


LniGon} “a piety ox ss apole fad 19480 


~ 


"{E} 


oon are! gist ork 2 


rik Ty Gee 9: ak a 
~sigeetoy: g3dq io. aokapest) git - erie j 
aa sini) net noone po + 74 , ee a : 
te rohass 12 iegew gest, 2 e 7 ee 
' SNe ae S705 a | Vaated C 
ek Te ee mn? 
to s2ettadah ne io ere. rs Y ognabvamanions = an , ® 
te: pofieokia ade, eae ae 5 Grout 


Dropping the capillary pressure term lead to a double valued 


water saturation at the oil-water interface which is a physical 
impossibility. To overcome this Buckley-Leverett propose a method of 
balancing areas at the flood front to establish a single sharp inter- 
face at the front. Controversy arose over this fact and Douglas 

eC. ait did a numerical integration of the Buckley-Leverett approach 
with the inclusion of the capillarity and gravity effects. Their 
results showed only a rounding off of the saturation distribution 


(Figure 3) at the flood front, the degree of rounding depending on the 


magnitude of the capillary and gravity forces. 


Rapoport and Dedeua elaborated on the Buckley-Leverett theory 
describing the transient flow phenomenon in a linear system. Rapoport 
and Leas assumed incompressible immiscible fluids in the absence of 
gravity forces and placing the expressions in dimensionless form, 


Rapoport and Leas solved five simultaneous equations to obtain: 


ds dF 0S K 1 ce) K. -F ‘dP0S 
Ww WwW ro Cc W 


— +——— - =| = 0 4 
oT TGs oX C_ LVy_ ox dS 0X (4) 
; Ww fe) Ww W 
when T = tV/Ld = dimensionless time co-ordinate 
Keo 5 
F= 1+ = dimensionless function 
C.K : 
o rw of saturation 
X = x/L = dimensionless space co-ordinate 
Ho 
C =-— = dimensionless oil-water 
4 By viscosity ratio 
V = total flow rate per unit cross- 


sectional area ( cm/sec ) 


and K = specific permeability ( darcies ) 
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for the boundary condition: 


K a7 1 +a : Kk, : dP os, = (0, at X = O for any T (5) 
a dS 9X 
Ww 

Equations (4) and (5) supplemented by a statement of the initial 
séturation distribution completely describes the flooding behaviour 
in a linear horizontal porous medium. This non-linear parabolic 
differential equation is not Only a function of the rock and fluid 
properties but also a function of the system length and injection 


rate as shown by the factor LVL, - However this factor does not 


account for the influence of capillary and wetting forces. 
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MODEL SCALING 
Before an expensive engineering project (such as a waterflood) 

is undertaken it is advisable to examine the performance of a small 

scale replica (model) of the reservoir (prototype) whose behaviour is 

to be studied. Thus it is necessary to define both experimentally 

and theoretically the extent to which the model is truly an analogy 

of the prototype. The design of a dimensionally sealed experiment or 

pure modek that will behave like the prototype involves, in part, the 

establishment of a number of dimensionless groups relating the various 


properties of the system. 


Rapoport and lease ue introduced a scaling coefficient LVu to 
study the effect on waterfloods of variations in length, velocity and 
viscosity. It was found that if the value of LVL was greater than 
one, then capillary end effects in the model would be negligible. 


ae that this is not a universal coeffi- 


However, it has been shown 
cient because it does not account for the capillary forces. Other 
similar groupings have been proposed in which parameters such as 

contact angle, porosity, interfacial tension and permeability have 


been introduced. 


A dimensionless scaling factor I was suggested by Engelberts 


and Klinkenbexg-°, deHaan’! and Goltine to scale the ratio of viscous 


to capillary forces, however, it neglects the contact angle and the 
porosity. 
LVu 
yee (6) 
y *K’ 
I = ratio of viscous to capillary forces 


L = total length of the flooded system (cm) 


Y = interfacial tension between the oil 


and water (dynes/cm) 


K = absolute permeability (om) 
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Although I is an useful scaling factor, Geertsma et ai showed the 


group to be a simplification of the universal rate function lL, 


LV ae 
Ar Fag ee menrpes eas (7) 
y cos68 v K¢ 
8 = contact angle made by the water with 


the oil on the rock surface 


This later group includes the variables of porosity and contact angle 


which had previously been neglected. 


20 ; F , 
Chuoke et al , carried out a theoretical and experimental 
analysis of fingering. A dimensionless mathematical group was de- 
veloped which related the average distance between fingers to various 


properties of the system. 
K 
Wen yee oe 
m CS Te, (8) 


da = the average distance between fingers(cm) 


The proportionality constant, C, is a dimensionless constant 
experimentally obtained and was found to have a value of 30 for a system 


with no connate water present. 


The relationship between a and the smallest finger thickness 


possible at the onset of fingering is given by: 


Bs feces ((S) 


Ay is called the critical fingering distance. The critical fingering 
distance as a fraction of the tube diameter in a linear flow system 


was found to be 


r = 1 
aa =C (14. (10) 
Say (uy = rae) h 


is as given is actually the inverse of what was proposed by Geertsma 
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h = the diameter of the core (cm) 


i Aah is greater than unity no fingering occurs, but if 
the value is less than one a large number of fingers develop. By 
solving equation (10) at the onset of fingering ( any =) 1) che 


value for C may be calculated. 


Vas’ Usa h’ (11) 
(6 = ——————————————EEE 
YK 


Since the value for C has been found to increase if connate water is 


present it must be established experimentally for each system. 


By solving equation (7) for y and substituting into equation 
(10) a relationship between L and the critical fingering distance 


may be obtained: 


ae es ae (12) 
m L. cosé ¥ o'( Uo / My -1) he 


and the value of L. may be calculated at the onset of fingering for 


each porous medium. 


DeHaan, Engelberts and Collins observed a decrease in break- 
through recovery with an increase in the scaling factor once the 
fingering commenced. By calculating the value of L for the point at 
which fingering occurs ( A ayfh = 14) Lom equation(10), it was found 
to correspond with the point of decreasing recovery with increasing 


values of the scaling function. 


2k Jod ,etyo50 gakiegnti on x ea 


yi .qoiswab Pci ‘Yo, 0, sedan 
Ee 
af1 ( l= a\, aise 


og 
* 
“ 


(£1) 


ai Tetaw eismboo ht sasezont 03 “Sauo? iu 


" mataye ni ge ve aa 


soltpups ose : eee 
g0nBIerb aniysant? TED 


wads 9280 tee 
ts Jntog ‘ds: 10% ra 49, 


' 
at. a's 28) ar 
ve = eal " My he 
‘ ol a wd va eu Wee 
cn i 1 Re any 
. ' 
ey - ; 4 
; - ; ’ 7 | 
y hs ao \ 
' © 1 wear 2 ce : A 
. iY G 
} . La P 
bade i 
7 : . =) 
& 
3 5 | i <2 . if 
t i ae re 
q re rer, 
8 I i : ; : o 7 
= t* ’ ; 


APPARATUS AND MATERIALS 


A schematic diagram of the experimental equipment used in this 
study is shown in Figure 4. The nucleus of the apparatus consisted of 
a three part core holder into which the test sand was packed. Each 
core section was of high strength lucite with an inner diameter of 
4.08 cm, a 0.635 cm wall thickness and a length of 134.78 cm. These 
segments met firmly at the mid-point of threaded couplings and were 
sealed by means of O-rings and high vacuum silica grease. The total 
length of the core could be varied by attaching or disconnecting sec- 
tions of the flow tube. The optical properties of the lucite ensured 


visual observation of the flood front. 


Each end of the core was fitted with a lucite cap or endplate 
(Figure 5) that had a recessed interior and a threaded wall inset with 
O-rings. The cap accommadated a spreader plate to disperse the fluid 


and a 200 mesh stainless steel screen to contain the sand. 


Two brass electrodes located directly behind the concentrating 
plate and wired to an electronic relay, timer and an electric bell 
were employed to detect water breakthrough at the outlet end of the 
core. Upon completion of the electrical circuit by the water, the 


bell and timer were triggered. 


Pressure taps were installed just beyond the endplates and con- 
nections made to a Validyue transducer, Whittaker carrier demodulator 
and an Advance digital multimeter recorder. High strength copper or 
stainless steel tubing at the inlet and outlet end carried liquid 
and gas to and from the core. The core assembly rested on a steel 
stand and the angle of dip for the core could be set by varying the 


height of the system at the downstream end. 
Imbibition tests were conducted with an apparatus similar to 


the one described by Bobek, Mattax and nenekas~ Surface and inter- 


facial tensions were measured with a Cenco DuNouy tensiometer; and 


15 


-AamaseiTA 


eidi «wt hoav siden (asin 
Io boteteron' auIstnage aia tn i 
dish Agateg, em) boas Jaad ody 

te sesouats senint ry pee 
snadt mo 4 Aeie 20 Sattol 8 ‘has e 
attow’ bre agnttqeo bobusads "aa Bes 
Lagat SHT Jopeaty vy a3tL Ds sine dati f 
~598 gaisosdnonets ‘te Srridoet ta) ad. obipe so bi 
boywwens Sat: as], ody is ookrgeqoss, [notsqo, . 
aoe 


oh sibs | 16 qéD oitaul « Asse ba33 
div Joaqhetiow bsbsotdl3 bas x 
biol? Sid. se teqeth oJ saald ak 


see Sis iieinos “Od 2 moetoe I 


o " , 


ene 


ada 46 ae 3553ue ala i one 
eit etetaw a3 “as sia fssinioat 


‘xeielgbomeb. rere ba eee 


“0 sealrondaige 3, agin “aeptoo0x Breese 
bitrip tf be} r19 bre aatsoe bas toink baad 18) sate 


be inbox zen ata a5 coels at. He8 
? 7 | ‘i 3 oe * 
os s8if Its GHIBTN a8 My ‘tt itw aovniion ogy ned), 


= + “youn! bog sveite . lane ben 8d 


Hires, i 5 St sia es yuoviod. ‘pond B 3 tw, bowense oT3t. anes takos 


af , 5 
‘ ; ‘ ; : Ps 


wa 


16 


SNLVEVddVY IWINSWINSdxX3S SHL 4O OILVWAHOS :% 3YNSIS 


dWNd eS 
vysny aivyu| 8 
YRGNINAD : LNVLSNOD OL] ~ 
3.1vNavUs 
300419373 ONITdNOD dv) ON3 
03378 
QOOH mi 
3WN4 OL eS oh 
YOLVINSaY S 
3unssaud & a > 
WOVE Ss Wind | 3 
On uouvinaowaa[___ 
Av ay VILLLL 
JINOY 19313 y¥30y4093u Alddns| _, 
YSLAWONVW-| lv | 3 
OL | 8 
1 E : 
38 x Yai z 
Y3NIVLNOD 
yaVM 


NOTIVS 3AlJ 


ine 


. 
b 
OB 
ms) 
> 
ro 
: 
ta 


: 
ae 
Le) 

uw 
-o aa 
vd 

a 
QO | 
=. 
Pal 


5 


é 


' 


d¥VO QGN3 
3H1 43O NOILYOd G3SS393Y¥ OLN! Lid N33YOS 
ONV 31V1d Y¥30V3S"dS SHL S30d0¥LI313 OML 
SHL HLIM dVD GN3 L37LNO 3HL 4O NOILO3S 
-SSOYD V ONV 31V1d Y30V3SYdS 3H1L *N33N0S 
7331S SSAINIVLS HSSW 002 SHL 40 SWVYSVIO :S 3YNSIS 


| 


WL 


N 


vi 


o 


rk 


LiOv 


c 


4. 


” 
we a 


Seko 
LHE CbBEVDED br Vie 


& 


‘ 


~ 
i] 


a 


EL WO HEC 
SECLION OL LRE 


“£MO  EFECL BODE 


“2CwEEM 


CVb MiLH 
v 


‘S 


WE 
SSL 


vi 


eM 


eae Py 
SHEVDEW Br VL 


E 


es 
wii 


| ae 
2eOt LAE 


= 


-2CMEEM’ LYE 


_TI@AwE 2 * -mvesy 


ba 


18 


the crude oil viscosities were evaluated from measurements obtained 


with the Bendix Lab Viscometer. 


De-aerated distilled water was chosen as the displacing fluid 
in the waterflood tests. The water was coloured yellow with uranine 
dye to permit visual observation of the displacement front. Uranine 
dye was selected because: 

1. No reaction takes place between the solute and 

the grains or the solute and the crude, and 

2. The solute imparts flourescent properties to 

the solvent which facilitates visual and optical 

observations of the flow behaviour. 
Iso-octane and Pembina crude were selected as the oils to be used in 
the displacement tests. Propane in conjunction with carbon dioxide 
was used to flush the oils from the core during the cleaning process. 
A Grove valve and back pressure regulator was attached at the outlet 
end to hold a back pressure so that propane could be kept in the 


liquid state. 


The,test sand was a clean, unconsolidated Ottawa sand 
(Fisher Scientific S-151) which was passed through an 80 mesh screen 
and retained by a 120 mesh screen. Examination under a binocular 
microscope disclosed the sand to be composed almost entirely of sub- 
rounded quartz grains with minute ferromagnesium impurities. A 
nearly constant grain size distribution was used to ensure reproduc- 


ability in repacking different cores. 
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PROCEDURE. 


To minimize any possible bypassing at the wall of the lucite, 
the inner surface of the core barrel was coated with a thin layer of 
sand. This was accomplished by filling che lucite holder with 
ethylene dichloride for approximately ten minutes, emptying the fluid 
and filling the tube with Ottawa silica sand. After approximately two 
minutes the sand was drained leaving a thin coating of quartz grains on 
the core barrel wall. It was observed that the lucite developed fine 
cracks after being exposed to the ethylene dichloride. An examination 
disclosed the crazing to be only surface deep and therefore would have 


no effect on the core columns strength. 


The sand was packed by agitating the core in a vertical position 
with high frequency vibrators while the grains fell through a two foot 
head of water at approximately one centimeter per minute. Vibrators 
were placed below, above and at the water-sand interface. Following 
packing, the core was vibrated and flooded simultaneously for a day to 
check for any sand settlement. The column was then dried by injecting 


air through the core for three days. 


Preceding water saturation of the core, a vacuum corresponding 
to a 27 inch drop in mercury was pulled on the system. Water was hydro- 
statically fed into the core until it began producing at the outlet end. 
The vacuum system was then disengaged and a constant rate Ruska pro- 
portioning pump was used to inject. water into the core until a steady 
state flow rate was established at the downstream end. A material 
balance calculation, based on careful measurements of the inlet and 
outlet volumes of water, yielded the porosity. Next, pressure drops 
were recorded at different flow rates and by solving Darcy's law the 
absolute permeability was calculated. The core fluid was then displaced 


by oil to an interstitial water saturation before a waterflood was run. 


At the outlet end a 2000 ce and twelve 100 cc graduated Kimax 
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cyclinders were assembled to collect the production of fluids. The 
larger cyclinder was used to collect oil until water breakthrough and 
the smaller cyclinders were them employed to collect fluids for the 
remainder of the waterflood. The 100 cc graduates were changed at 
specific times and the volume of water and oil in each was recorded. 


Using the data, the water-oil ratios were calculated. 


Before an oil was displaced by water the core was cleaned in 
the following manner: 
1. The sand pack was flushed with liquid propane 
until no further oil was produced. 
2. co, was injected to displace the propane and 
partially dry the core. 
3. Air, fed through a dessicator, was then in- 
jected through the core for 2 to 3 days to 
completely dry it. 


Upon the completion of the run #12 one third of the core was 
removed and experiments conducted on the remainder of the model to 
check the effects of length on recovery. The segment of the column 
removed was repacked with the original sand to examine the reproduc- 
ability of the initial wetting and recovery conditions in the porous 


medium. 


Imbibition tests were carried out on mini-cores with a length 
of 3.81 cm and a diameter of 1.53 cm. The cores were placed in a 
vacuum dessicator and saturated with the desired fluid. Once saturated, 
the core was placed in the imbibition apparatus, surrounded with the 
imbibing liquid and left for 24 to 72 hours during which time the fluid 
imbibed into the sample. Imbibition tests were completed: 

1. the original sand saturated with water, 

2. the original sand saturated with oil, 

3. an oil saturated sand taken from the 


core after run #12, 
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4. the sand, in the core after run #12, which 
had been cleaned, dried and resaturated with 


oil, and 


5, the sand, in the core after run #12, which 
had been cleaned dried and resaturated with 


water. 
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DISCUSSION OF RESULTS 


Values of the initial rock and fluid properties (¢, K, u, Y ) 
before a waterflood, were measured by means of preliminary tests (tests 
conducted prior to oil displacement). Table 1 gives a tabulation of 


these findings and shall be expanded upon below. 


The methods for the evaluation of these parameters have been 
previously discussed in the section titled "Procedure". The values for 
the porosities varied from 38.7 to 39.7 percent displaying a deviation 
from the mean of + 0.50 units. This is considered to be an acceptable 


result well within the range of experimental error. 


Permeability was calculated using the classical approach for the 
evaluation of Darcy's Law. Between run #5 and #6 the core was cleaned 
and a change in the absolute permeability from 16.8 to 15.2 darcies 
occurred. This drop was attributed to the incomplete cleaning of the 
core and/or possible incomplete resaturation with water between runs. 
Packs for run #15 and #16 yielded slightly lower permeabilities of 15.2 
and 14.5 darcies respectively. For these much shorter cores it was 
postulated that the closer spacing of the three high frequency vibrators 


lead to a slightly more consolidated pack. 


Iso-octane and Pembina crude were found to have viscosities of 
0.58 and 5.96 cp respectively at room temperature. Interfacial tensions 
between each fluid and water at room temperature were 41.55 dynes/cm for 
iso-octane and 32.92 dynes/cm for Pembina crude. These values were 


considered to be normal for the above oils. 


Table 2 shows a tabulation of the core and fluid characteristics 
present in the sand pack before and following each waterflood. Iso- 
octane ( Woy, < 1) was the displaced fluid for tests #1 to #5 while 
the remainder of the tests employed Pembina crude ( Hof ee, Pj. Initial 
water saturations (S.4) for the iso-octane floods ranged from 26 to 28 


percent whereas for the Pembina crude floods initial water saturations 
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varied from 20 to 27 percent. A decrease in the initial water 
saturation after changing from a lower viscosity to a higher 


viscosity oil was also observed by Coudinees 


To test the reproducibility of results, tests #1 and #2 
were conducted under similar conditions and yielded breakthrough 
recoveries of 72 and 71 percent; at a WOR of 6 the oil recoveries 
were equal at 73 percent. With these findings it was concluded 
that the results are within experimental error and indicate an 
acceptable reproducibility of a test. Further information can 
be obtained on the recovery histories for each run in Appendices 
B and C. Appendix B gives a tabulation of this data and 
Appendix C includes a graphical illustration of recovery versus 


pore volume of water injected. 


Figure 6 illustrates the recovery versus pore volume of 
water injected for typical cases of the two oils employed. 
Referring to the graphical representation of the recovery history 
for iso-octane, it may be observed that the production obtained 
up to breakthrough is 72 percent. After breakthrough subordin- 
ate production is almost zero. This result closely agrees with 
the work of Collins and Rapport et he for similar viscosity 
oils. The recovery breakthrough constitutes the major fraction 
of the ultimate oil eecoreey due to the favourable mobility 


ratio of the system. 


The recovery curve for Pembina crude in Figure 6 
illustrates a lower breakthrough recovery of 44 percent but a 
more significant subordinate production. This is again due 
to the mobility ratio for the system. The larger viscosity of 
the Pembina crude leads to a less favourable mobility ratio 


resulting in a lower breakthrough recovery and greater sub- 


ordinate production. 
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From the tests it was concluded that the wetting conditions 
changed in the sand pack while Pembina crude was being used in the 


core. Three distinct phases were postulated: 


1) a water-wet case, 
ii) a transition or possibly neutral case and 


7 af hp an oil-wet case. 


Justifications for making this deduction are based on the following 


evidence. 


The values of effective permeability to oil (K,) as tabulated 


in Table 2 may be divided into three categories. 


i) K, values of 5.0 darcies (runs #6, #7, #8, 
#16, and #17) were interpreted as the water- 


wet case, 


sky 2 Ky values of 3.0 darcies (runs #9, and #14) 


were interpreted as the transition case and 


a ba is By Ko values of 2.0 darcies (run #10, #11, #12 


and #13) were interpreted as the oil-wet case. 


The Welge Method as outlined in Appendix D was used to 
calculate the permeability ratios from measured production data. 
Figure 7 represents the relative permeability curves calculated in 
the above three cases. The shift in the effective permeability 
ratio from A to B to C is characteristic of changing wetting 
conditions. However, one anomalous curve was obtained from run #7. 
All tests conducted with Pembina Crude had an initial water satur- 
ation between 20 to 24 percent except for run #7 which had an 
initial saturation of 27 percent. Fatt and giikierece demonstrated 
that a shift of the relative permeability curve to the left will 
occur if the porous medium changes from a water-wet to an oil- 


25.45 : 
wet system. However, Stewart in his research found that 
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TEST SYMBOL 


WATER-WET (Swi =27 %) 


‘20 30 40 50 60 70 
WATER SATURATION ( PERCENT) 


FIGURE 7: RELATIVE PERMEABILITY CURVES FOR 
THE THREE WETTING CONDITIONS DISPLAYED 
BY PEMBINA CRUDE. 
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an increased initial water saturation shifts the relative permeability 
curve to the right. It was therefore concluded that the anomaly 
shown for the curve representing run #7 was due to the much higher 


initial water saturation. 


Imbibition tests were conducted on a sample of the original 
sand and a sample of the sand taken from the core following run #12. 
A complete description of the results for each test, is given in 
Appendix A. From these tests it was concluded that the original sand 
was water-wet and that the sand sample taken from the core after pro- 


longed contact with Pembina crude was oil-wet. 


With the exception of run #8 the average breakthrough recovery 
for the tests conducted where the effective permeability to oil values 
were around 5.0 darcies, was observed to be about 43%. For runs #10 
to #12, where the effective permeability was about 2.0 darcies, the 
breakthrough recovery was found to be about 22 percent. This trend 
was also observed by deuaan’ who had controlled wetting conditions, 
and is characteristic of the water-wet and the oil-wet cases 


. 


respectively. 


The previous explanations for assuming a change in wettability 
are by themselves not conclusive. However, when these reasons are 
taken in conjunction with one another they form a reasonable basis to 


substantiate the wettability change. 


Table 3 gives a tabular representation and Figure 8 is a 
graphical illustration of the calculated values of I as a function of 


recovery for the various tests, where I may be found as follows: 


Evi 
I =—— (6) 
y Vik? 
and Li: 
LVuU 
if be (7) 
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TABLE 3 


SCALING FACTOR VALUES 


(1) 


Breakthrough 
Recovery 


43 


39 


30 
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An examination of Figure 8 shows the recovery of oil at breakthrough 
(during the waterflooding of iso-octane) increase as the scaling 
factor increases to a value of about 0.5. Above this value the 
recovery is fairly constant. The observed decrease in recovery 
below 0.5 is attributed to the capillary force end effects which play 


a dominant role at low displacement rates. 


The waterfloods conducted on the Pembina crude gave data which 
resulted in the three breakthrough recovery versus scaling factor rela- 
tionships shown in Figure 8. For the water-wet case a noticeable 
decrease in recovery was observed when I exceeded a value of 
approximately 3.0. Below this value of I, recovery was stable to a 
value of about 0.5 where a decrease in recovery occurred due to 
capillary forces. Above a scaling factor value of 3.0 the decrease 
in recovery is attributed to the onset of viscous fingering of water 
into the oil which results in isolated oil pockets being left behind 


the flood front. 


The breakthrough recovery for the transition and oil-wet cases, 
of 30 percent and 23 percent respectively, in Figure 8 was found to 
be fairly stable over the range of rates studied. However, it is 
possible that viscous fingering could occur under these wetting 


conditions for very large throughput rates. 


In test #13 the core length was decreased and the flood carried 
out at the same rate as run #11 -(614 cc/hr). The breakthrough 


recovery in both tests was found to be 23 percent. 


When employing the scaling group I it must be remembered that 


the group is not universal but only a specific form of the Ly function. 


: 3 27 
Figure 9 shows a graphical representation of the studies of Engelberts 


and deHaan who related the breakthrough recovery to the same scaling 
28 

group I. By considering the previously mentioned work of Chuoke  , 

the onset of fingering was mathematically determined for the work of 


Engelbert and deHaan. The calculated value of I for ,,=h (the point 
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below which fingering will occur) corresponds to the point of decreas-— 


ing recovery with increasing I values as shown in Figure 9. 


Both Engelberts and deHaan used similar unconsolidated cores 
in their experimental apparatus. These cores were originally saturated 
with oil and would only imbibe water to a slight degree. This was 
referred to by deHaan as an 'intermediate-wet ' system. Since their 
cores were not truly water wet, the contact angle should not be 
neglected (as deHaan and Engelbert did in evaluating the scaling factor), 
otherwise a comparison can only be made with another intermediate-wet 
porous medium. Justification for omitting the porosity, as was done by 
Engelberts and deHaan, is only valid if the percent pore space in the 
two systems to be compared is equal. Since the systems of deHaan and 
Engelberts were nearly identical and experimental conditions were very 


similar, this may be the reason for their close agreement. 


A comparison of Newcombe tea: and deHaan's work is shown in 
Figure 10 and a very close agreement exists. This agreement may be 
attributed to, the fact that the contact angle was the same for both 
cases. Newcombe's experimental procedure for flooding silica sand was 
the same as that of Engelbert and deHaan. In his work Newcombe showed 


that the contact angle at breakthrough was 83°. 


For the purpose of comparison with this work, a contact angle 
of 83° was assumed for Engelbert's and deHaan's work. It is necessary 
to assume a suitable contact angle for this study and Collins' work so 
that all data may be compared using the universal scaling factor Li. 
teseho” reported a contact angle of 30° for Berea sandstone for tests 
conducted under conditions similar to this study and the work of 


Collins. Therefore, 30° was used as an approximate value for the 


water-wet case for these two studies. 


The values of L. for the work of deHaan and Engelberts was cal- 
culated assuming a contact angle of 83° and corresponding porosities of 


33 and 38 percent respectively. For the work of Collins and this study a 
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contact angle of 30 and porosities of 36 and 39 percent were used. 
Figure 11 illustrates the scaling factor L. versus breakthrough re- 
covery curves for the above four cases and from these results the 


following observations were made: 


i) For the similar viscosity ratios of this study 
and the work of deHaan the point where decreas- 
ing recovery with increasing L. values occurs 
and the calculated values for the initiation 
of fingering ( ee = h) show close agree- 


ment to one another. 


ii) There is a slight shift of the results to the 
left of oe = h for the larger viscosity 


ratios of Engelberts and Collins. 


iii) In all cases the onset of viscous fingering 
appears to occur at values of L. greater than 


1. Os 


It is therefore, safe to assume that no fingering will occur if L. 
is less than one in an isotropic homogeneous reservoir as long as 


the viscosity ratio is not greater than that used by Collins. 


The variation of recovery for the four experimental curves 
shown in Figure 11 is a function of the difference in mobility 
ratios and the pore geometry for each system. However, due to 
insufficient information the data is represented in the illus- 
tration as a function of the fluid viscosity ratio rather than the 


mobility ratio. 


If a value for the proportionality constant C is known, a 
value of the scaling factor L. for the onset of fingering may be 


evaluated from the following equation: 
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For a porous medium containing no interstitial water C has been found 


to be equal to 30. 


In this work the value of C was calculated to be 64. Hence, 
for isotropic homogeneous reservoirs, where the value of L, K, 
contact angle, Uo» Woosh, > and Y can be estimated, an approximate 
value for velocity of the flood front, beyond which viscous finger- 


ing will occur, may be established. 


By assuming a linear homogeneous isotropic water-wet porous 
media having properties similar to the values in the table below, it 
was found that under reservoir conditions the velocity (V), of the 
advancing waterfront, must be less than 4.2 x ne ft/day if no 


fingering is to occur. 


Length 300 meters 
Water Viscosity ep 
Interfacial Tension 30 dynes/cm 
Contact Angle 30 degrees 
Permeability 100 md 
Porosity 30 percent 


Once a value of 4.2 x 191s ft/day for the flood front advance is 
exceeded the chances of viscous bypassing are greatly increased. 
This implies that at least within the vicinity of the injection 


wellbore fingering may occur due to the high velocities. 


In the laboratory it is much more common for the finger 
height to be greater than the height of the core model. In this 
case a stable frontal displacement occurs in the model leading to 


higher recoveries than in the reservoir. 
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CONCLUSIONS 


Based upon this experimental investigation the following 


conclusions have been reached: 


a 


For an unconsolidated sand pack exhibiting water-wet conditions 
the stabilized breakthrough recovery as a percent of the oil in 
place was 43 percent and 72 percent for viscosity ratios of 6.4 


(Pembina crude) and 0.6 (iso-octane) respectively. 


The adsorption, over a prolonged period of time, of what 
appeared to be a surface active constituent from the 
Pembina crude onto unconsolidated quartz grains resulted 
in a wettability change of the Ottawa silica sand from 


a water-wet to an oil-wet system. 


For values of the scaling factor Le less than one, in a 
water-wet system, a stable front forms and recovery is 
essentially independent of a variation in the scaling 


factor. 


The value of the proportionality constant C was found to 
be 64 for the unconsolidated Ottawa silica pack used in 


this study. 


Although there are other factors that might effect viscous 
fingering over the range of rates studied, no fingering 
occurred when the porous medium was in a neutral or an 


oil-wet state. 
When the mobility ratio was less than one (iso-octane) no 


viscous fingering was observed over the range of the rates 


studied in this investigation. 
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-RECOMMENDATIONS | 


To further study the effects of the physical parameters of the 


scaling factor Ly on oil recovery the following tests are recommended: 


i. Use of unconsolidated cores with varying grain size 
distribution to act as a check on the effects of 


porosity and permeability and 


ii. Conduction of tests on consolidated cores, both 
homogeneous and heterogenous in nature, to further 


evaluate the universalism of this scaling group. 


With the aid of refined oils of varying viscosities it would 
be possible to obtain a better understanding of the change in the 
onset of fingering ( do = h) with oil viscosity. Refined oils would 
tend to eliminate the wettability change that occurs when using natural 


crude. 


In order to properly evaluate the effect of initial water 
saturation on the proportionality constant C it is recommended that 


the initial water saturation be varied on the same test core. 


Measurements of contact angles on different rocks for varying 
oil viscosities would prove invaluable. This would allow future 
researchers access to a contact angle representative of their system 


when calculating the onset of viscous fingering. 


By flooding the system at different angles of dip the recovery 
performance and the critical velocity for the onset of water tongues 
as predicted by Dietz could be compared to the actual measurements. 
At higher oil viscosities (HU, 7 10) the use of the viscous fingering 
method of van der Poel and van Meurs versus the Buckley-Leverett 
theory for predicting performance behaviour could be compared in the 


region where viscous forces play a dominant role. 
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LIQUID AND CORE PROPERTIES 


In order to establish the surface and interfacial tension of 
water, Iso-octane and Pembina Crude ten readings were taken on the 
duNouy tensiometer at room temperature. The arithmetic average was 
then calculated and corrections applied to thé mean which yielded the 


following results: 


SURFACE TENSION INTERFACIAL TENSION 
BETWEEN OIL AND WATER 
(dynes/cm) (dynes/cm) 
ISO-OCTANE 19.95 AE We p>) 
PEMBINA CRUDE 2Ia33 32.92 
WATER TREO? 


Viscosities of Iso-octane and Pembina Crude measured with the 
Bendix Lab Viscometer at room temperature gave results of 0.58 cp and 
5.96 cp respectively. The viscosity of distilled water at various 
temperatures was obtained from Perry 'Handbook of Chemical Engineering" 


and Figure A-1l was constructed. 


The core consisted of three pieces of high strength lucite 
134.78 cm long and 4.08 cm in diameter. For runs #1 to #12 the core 
length was 404.34 cm long. During runs #13 and #14 the length was 


reduced to 269.56 cm and for run #16 and #17 the core was 87 cm long. 


The porosity was obtained by applying a material balance across 
the core while the initial saturation was being established. Reproduce- 
ability of porosity was satisfactory as shown by the following values 


which were established after each new sand pack was prepared. 
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Imbibition tests conducted on the sand gave the following 


information: 


DESCRIPTION 
OF SAND 


-original sand 
saturated with 
water 


-original sand 
saturated with 
oil 


-oil saturated 
sand taken from 
the core after 
test #12 


-sand taken from 

the core after 

test #12 which 

had been cleaned 

and saturated with oil 


-sand taken from 
the core after 
test #12 which 

had been cleaned 
and saturated with 
water 


TIME IN IMBIBITION 
CELL 


24 hours 


24 hours 


48 hours 


72 hours 


72 hours 


COMMENT 


-no water produced and 
no imbibition of oil 
into the core 


-some oil produced and 
water was found imbibed 
in the core 


-no oil was produced 
and no water was found 
in the core 


-no oil was produced 
and no water was found 
to have imbibed into 
the core 


-some water was displaced 
and some oil imbibed 
into the core 


This lead to the conclusion that the sand pack was originally 


in a water-wet state but due to prolonged exposure to the Pembina 


Crude changed to an oil~wet condition. 
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DISPLACEMENT TEST RECOVERY DATA 
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DETERMINATION OF PERMEABILITY RATIOS 
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DETERMINATION OF PERMEABILITY RATIOS WELGE METHOD) 


Employ the Buckley-Leverett equation under the following 
assumptions (approach is that proposed by Collins, R.E.) 

1) capillary and gravity effects are negligible 

2) fluids are immiscible and incompressible 


3) the porous media is homogeneous and linear 


For a system of length L with s (x,0)=S introduce a wetting 
fluid which displaces a non-wetting fluid at a rate q(t) at x=0 for 
t>0. At x=0 only the wetting fluid is flowing and at x=L both fluids 
are flowing. At the outflow end the cumulative production of the non- 
wetting fluid is designated by Oe and the inflow of wetting fluid is 


shown by Q, where 


crt a Be 


At the outflow end 


-1 
dQ vi 
y 2 Gow _ nw _ 1 + bw D - (2) 
nw q dQ K wen 
or from the definition of Ee 
K u res sl ast 
7 ow Dw. to nw Jad D—- (3) 
K U dQ dQ 


Welge developed an equation for calculating the non-wetting 


saturation: 


periogan (i supra. Seregde 
S g(t) “all a eee x D- (4) 


dQ 


A graph of Qaw Versus 0: 
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Ow 
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shows that from Q = O to Q, the line has a slope of 1 representing 


a production of oil only. For Q > Qe» then 


Qe. =a+tb 1nQ De 6) 


for any small segment. By differentiating we find: 
dQ 


b D - (6) 
for each segment. This may be used in equation D-(3) and D-(4) 
along with the value of OF to compute the permeability ratio and 
its corresponding S wil): A sample calculation is shown for run 
#1 in Table D-1 and graphed in Figure D-1. Points from the 
relative permeability curve were applied to equation 2 (for the 
conditions of oP /du = 0 and sina =0) to solve for the fractional 


flow of water at different water saturations (Figure D-2). 


From the fractional flow curve a theoretical water break- 


through time (T,) 


r AL} (Sa -Ss.,) 


wi 
Ce eee D- (7) 
Q 


was calculated and the results are tabulated in Table D-2. The 
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TABLE D-1 


SAMPLE RELATIVE PERMEABILITY CALCULATION 


RUN #1 


From the recovery graph of Oil Produced (Q) and Water Injected 
(Qaw the following values of Kank were obtained: 


Q dQ. She S ee 
dQ a 
1720 [555 495 505 76 
1750 "= 4233 1325 675 185 
1800 0667 234 P66 044 
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theoretical and observed times for rates of 120 to 1200 cc/hr yielded 
a maximum deviation of 5.8 percent. However, in the case of run #4 
(30 cc/hr) the error increased to 9 percent. This is due to the fact 


that capillary forces play a dominant role at low rates. 
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